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ABSTRACT 

The central engine of gamma-ray bursts (GRBs) is believed to be a hot and dense disk with hyper- 
accretion onto a few solar-mass black hole. We investigate where the magnetorotational instability 
(MRI) actively operates in the hyperaccretion disk, which can cause angular momentum transport 
in the disk. The inner region of hyperaccretion disks can be neutrino opaque, and the energy- and 
momentum-transport by neutrinos could affect the growth of the MRI significantly. Assuming rea- 
sonable disk models and a weak magnetic field B < 10 14 G, it is found that the MRI is strongly 
suppressed by the neutrino viscosity in the inner region of hyperaccretion disks. On the other hand, 
the MRI can drive active MHD turbulence in the outer neutrino-transparent region regardless of the 
field strength. This suggests that the baryonic matter is accumulated into the inner dead zone where 
the MRI grows inactively and the angular momentum transport is inefficient. When the dead zone 
gains a large amount of mass and becomes gravitationally unstable, intense mass accretion onto the 
central black hole would occur episodically through the gravitational torque. This process can be a 
physical mechanism of the short-term variability in the prompt emission of GRBs. Finally, the origin 
of flaring activities in the X-ray afterglow is predicted in the context of our episodic accretion scenario. 
Subject headings: accretion, accretion disks — black hole physics — gamma rays : bursts — instabil- 
ities : magnetic fields 



1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic 
event in the universe. GRBs are generally considered 
to be powered by hyperaccretion of the matter onto a 
stellar-mass black hole (~ 2-3M Q ), which is formed in 
the context of the "collapsar" scenario or merging sce- 
narios of compact objects (Narayan et al. 1992; Woosley 
1993; Paczynski 1998; MacFadyen & Woosley 1999; Ruf- 
fert & Janka 1999, 2001; Meszaros 2002; Piran 2005). 
The hyperaccretion rate is of the order of 0.1-1 MqS~ x 
and the release of gravitational energy powers the burst. 
The radiative energy ejected through relativistic jets is 
expected to account for the observed 7-ray emission. 

The key process for releasing the gravitational en- 
ergy is angular momentum transport in the disk (Lee 
& Ramirez- Ruiz 2006). As in the cases of the other 
astrophysical accretion disks, magnetic turbulence initi- 
ated and sustained by the magnetorotational instability 
(MRI) is believed to play an essential role for the angu- 
lar momentum transport in hyperaccretion disks (Balbus 
& Hawley 1991; Lee et al. 2005; Fujimoto et al. 2006; 
Nagataki et al. 2006). Detailed linear and nonlinear 
analysis of the MRI facilitate our understanding the role 
of the MRI in various accretion systems (Hawley et al. 
1995; Machida et al. 2000; Turner et al. 2002; De Vil- 
liers et al. 2003; Sano et al. 2004). However, the physi- 
cal conditions of hyperaccretion disks are quite different 
from the other accretion systems. Because the hyper- 
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accretion disks are very dense and hot like supernova 
cores, they can be cooled through the neutrino radiation 
(Popham et al. 1999; Narayan et al. 2001; Di Matteo et 
al. 2002). In addition, the energy and momentum are 
mainly transported by neutrinos in the neutrino-opaque 
region. However, the detailed analysis of the MRI in hy- 
peraccretion disks including such microscopic effects has 
not been done yet except for Araya-Gochez & Vishniac 
(2004). 

Masada et al. (2007; hereafter MSS07) have investi- 
gated the effects of the neutrino transport on the MRI 
in proto-neutron stars (PNSs), and shown that the heat, 
chemical, and viscous diffusions caused by the neutrino 
transport have great impacts on the MRI. The heat and 
chemical diffusions can reduce the effect of stratifications, 
and the neutrino viscosity suppresses the growth of the 
MRI. Because the resistivity of hot nuclear matter is rel- 
atively low (Yakovlev & Shalybkov 1991; Miralles et al. 
2002), the magnetic diffusion is much smaller than the 
other diffusion processes, and hardly affects the MRI in- 
side PNSs. 

The hyperaccretion disk is considered to have similar 
physical properties to the PNS. In particular, neutrinos 
can be trapped in the inner dense region of hyeraccre- 
tion disks, and the energy and momentum are mainly 
transported by the neutrinos. Therefore, we can apply 
the results of MSS07 to hyperaccretion disks with little 
changes. 

The growth time of the MRI in the absence of vis- 
cosity is given by X/va, where A is the wavelength of a 
perturbation and va = B / '(4tt p) 1 ^ 2 is the Alfven speed. 
The growth of the MRI is suppressed dramatically if the 
growth time is longer than the viscous damping time ~ 
A 2 /v. where v is the kinematic viscosity (MSS07). Then, 
a large enough viscosity can reduce the linear growth 
of the MRI. Clearly, as the wavelength becomes longer, 
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larger viscosities are required to suppress the MRI. Be- 
cause the typical wavelength of the MRI is A ~ va/Q, 
the condition for the linear growth of the MRI can be 
written as 

v v\l 

where Re is the Reynolds number, and £1 is the angu- 
lar velocity. Here we choose va/Q as the typical length 
scale L and va as the typical velocity U. The Reynolds 
number is a good indicator for rapid growth of the MRI 
in hyperaccretion disks. In the dense neutrino-opaque 
matter, the kinematic viscosity via the neutrino trans- 
port would be quite large, so that the condition given 
by equation (fl]) may not be satisfied. Therefore, in this 
paper, we investigate where the MRI operates in hyper- 
accretion disks focusing on the neutrino viscous effect. 

The plan of this paper is as follows: In § 2, we construct 
reasonable disk models, and derive simple prescriptions 
for the neutrino opacity and viscosity. Linear growth of 
the MRI in hyperaccretion disks is studied in § 3 using 
the dispersion equation derived in MSS07. Based on the 
linear stability analysis of the MRI, we propose a new 
evolutionary scenario of collapsar disks that can explain 
the short-term variability in long GRBs in § 4. In § 5, 
we discuss a possible origin of flaring activities in the X- 
ray afterglow, which have been discovered by the Swift 
satellite. Finally, we summarize our main findings in § 6. 

2. THE PHYSICAL MODELS 
2.1. Neutrino Depth and Viscosity 

We need to estimate the amplitude of the kinematic 
viscosity in hyperaccretion disks to investigate where the 
MRI is active. The momentum transport by neutrinos 
plays an important role in the neutrino-opaque regions. 
In contrast, its influences can be negligible when the ac- 
creting gas is transparent to the neutrino. We treat, for 
simplicity, the neutrino transport with the diffusion ap- 
proximation. The kinematic viscosity is assumed to be 
zero in the neutrino-transparent regions. Note that the 
diffusion approximation fails when the neutrino depth 
reaches to unity. However, we focus mainly on the qual- 
itative features of the MRI in hyperaccretion disks, so 
that we adopt this simple treatment in this paper. 

First we derive the neutrino depth and kinematic vis- 
cosity via the neutrino transport. Both absorption and 
scattering processes of the neutrino contribute to the 
neutrino depth. Absorption onto the nucleons by the 
inverse of the URCA process gives rise to a large por- 
tion of the absorptive neutrino depth. The other absorp- 
tive processes, such as the inverse of pair annihilation, 
bremsstrahlung, and plasmon decay would be negligi- 
ble (Kohri & Mineshige 2002). The absorptive neutrino 
depth is, thus, described by r a b s — 4.5 x 10~ 39 T 2 X nuc pH , 
where T is the temperature, p is the density, H is the 
scale height of the disk, and A nuc is the mass fraction 
of free-nucleons (Di Matteo et al. 2002) . The scattering 
neutrino depth through elastic scattering of background 
nucleons is given by r sc = 3pn sc H ~ 8.1 X 10~ 39 T 2 pH , 
where the factor 3 assumes equivalent contribution for all 
neutrino species and k sc is the mass scattering coefficient 
(Tubbs & Schramm 1975; Shapiro & Teukolsky 1983). 

Total neutrino depth T to t is given by the sum of the 



absorptive and scattering neutrino depth, 

Ttot = T abs + T SC = 1.3 X IQ^T 2 pH . (2) 

Here X nuc can be assumed to be unity in the inner dense 
region of hyperaccretion disks. The diffusion approxima- 
tion for the neutrino transport can be used in the re- 
gion where the condition r to t > 2/3 is satisfied. In such 
regions, the averaged mean-free-path of the neutrino is 
defined by (A) = iJ/r to t, and the neutrino viscosity v is 
described by 

, = i^(A)«5.2xl0 12 rV 2 cm 2 sec" 1 , (3) 
15 pc 

where U v — 4ergT 4 /c is the neutrino energy density (van 
den Horn & van Weert 1984; Burrows & Lattimer 1988). 
Here <tb and c are the Stefan-Boltzmann constant and 
the speed of light, respectively. Spatial distributions of 
the neutrino depth and viscosity in hyperaccretion disks 
can be obtained from equations © and in so far as 
the disk structure is determined. 

2.2. Structure of Hyperaccretion Disks 

We adopt basic equations based on the Newtonian 
dynamics and assume a quasi-steady structure of hy- 
peraccretion disks. Various steady disk models have 
been proposed as the central engine of GRBs taking into 
account detailed microphysics and/or general relativity 
(Popham et al. 1999; Narayan et al. 2001; Kohri & 
Mineshige 2002; Di Matteo et al. 2002; Kohri et al. 
2005; Chen & Beloborodov 2007; Kawanaka & Mineshige 
2007). These models are constructed in the framework 
of the a-prescription of turbulent viscosity (Shakura & 
Sunyaev 1973). However, we now have an interest in the 
features of the MRI that causes the turbulent viscosity. It 
is, thus, necessary to construct disk models independent 
of the a-parameter. For simplicity, we adopt power-law 
models for the radial distributions of all physical quanti- 
ties. 

The disk surface density S(r) is one of the most im- 
portant quantities in constructing the disk structure. We 
assume a power-law distribution with an index q, 

£(r) = S r- 9 , (4) 

where f = r/r s is the distance from the central black 
hole normalized by the Schwarzschild radius. The 
Schwarzschild radius is given by r s — 2GMsn/c 2 — 
8.9 x 10 5 M 3 cm, where M 3 = M B h/(3M ) is a mass 
of black hole Mbh normalized by 3 M & . Here So is a 
reference value of the surface density at r = r s . When the 
baryonic matter of its mass ~ 0.1-1 M is distributed 
within ~ 10 7-8 cm from the central black hole, the av- 
eraged surface density is in a range of ~ 10 17-18 g cm -2 
(Woosley 1993). Thus, we choose the reference value as 
So = 1.0 x 10 18 /i; g cm -2 , where /s is an arbitrary pa- 
rameter. 

Another key quantity is the disk temperature, which is 
assumed to retain a power-law distribution with an index 
P, 

T(r) = T f- p , (5) 

where To is a reference value at r = r s . A much less 
extreme assumption would be that the thermal energy 
is ~ 1 - 10 % of the gravitational energy of the disk 
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(~ GMbh/t). Then the typical disk temperature is esti- 
mated as ~ lO 11-12 K at r = r s . The reference value is 
chosen as To — 4.3 x 10 n /T K, where fx is another arbi- 
trary parameter. This temperature is equivalent to the 
Fermi energy of the partially degenerate electrons (Lee 
et al. 2005). 

Assuming the disk in which the gas pressure dominates 
the degenerate and radiation pressure, the sound speed 
is given by 

Cs = (^pj 1/2 = 6.1 x 10 9 /t /2 ^" p/2 cm sec" 1 , (6) 

where Ub is the Boltzmann constant and m p is the pro- 
ton mass. In a gravitationally stable disk, the vertical 
component of the gravity in the disk is contributed by 
the central black hole. The hydrostatic equilibrium in 
the vertical direction determines the scale height of the 
disk, 

H = ± = 2.5 x 10 5 /T /2 ^3r- (p - 3)/2 cm, (7) 

where the disk is assumed to rotate with the Keplerian 
angular velocity, f2 = £l K — 2.4 x 10 4 M^ 1 f -3 / 2 sec -1 . 
Then the density structure can be evaluated as 

p = A = 2.0 x l^M T 1/2 M^r^-^ g cm- 3 . 

Ltl 

(8) 

Using equations (JSJ), (J7J), and ©, the radial profiles of 
the neutrino depth and viscosity are given by 

7-tot = 1.3 X 10 3 / E /2^2 P - 9 ) (Q) 

and 

v = 8.1 x 10 10 / s 2 /t^3"^ 3(1 " p)+2<? cm 2 sec" 1 . (10) 

The strength of the magnetic field is also an impor- 
tant quantity to determine the growth of the MRI in 
hyperaccretion disks, although it is highly uncertain in 
the context of GRBs. We assume that the pre-collapse 
core of a massive star has a uniform magnetic field with 
the same amplitude as strongly magnetized white dwarfs, 
that is B p w 10 9 G and has an average density of 
Pcorc ~2x 10 9 g cm -3 . Considering that the magnetic 
flux is conserved in the core-collapse phase (B oc p 2 ^ 3 ), 
the radial structure of the poloidal magnetic field in hy- 
peraccretion disks is obtained as 

B p , disk = 1.0 x lQ u f^f-^f B Mf<*&-*-*>'* G , 

(11) 

where fs is the arbitrary magnetic parameter. The 
Alfven speed in the disk is then given by 

v A = 2.0xl0 4 /s /6 / T 1/12 /BM 3 " 1/6 r (p - 29 - 3)/12 cm sec" 1 . 

(12) 

The power-law indexes of the surface density and tem- 
perature are determined by the thermal equilibrium in 
the disk. There are two main cooling processes in hy- 
peraccretion disks, which are the advection cooling and 
neutrino one. Assuming that the advection cooling dom- 
inates the neutrino one, the power-law indexes of the sur- 
face density and temperature are given by p = 1.0 and 
q = 0.5 (Di Matteo et al. 2002; Chen & Beloborodov 
2007). This corresponds to the ADAF (advection domi- 
nated accretion flow) type disk. We can also choose the 



power-law indexes as p = 1.5 and q = 0.0, if the neu- 
trino cooling is dominant. This is so called the NDAF 
(neutrino-cooling dominated accretion flow) type disk, 
but the neutrino is assumed to be thermal here (Woosley 
1993). When the inner dense region is opaque to the neu- 
trino, the disk structure would approach either of these 
two models. On the other hand, the disk structure may 
be quite different from them in the neutrino-transparent 
region (Popham et al. 1999; Narayan et al. 2001). How- 
ever, since our main interest is in the neutrino-opaque 
region, the disk models with single power-law distribu- 
tions are adopted in this paper. 

The linear stability analysis shows that a Keplerian 
disk is gravitationally stable only when Toomre's Q = 
Qcs/^ttGE) is larger than about unity (Toomre 1964; 
Goldreich & Lynden-Bell 1965). The disk model given 
by equations ^ - © is gravitationally stable when the 
inequality 

Q = 7.0 x 10 2 / s Vt 2jW 3 _1 ^" (p_2,+3)/2 £ 1 . ( 13 ) 

is fulfilled. Since we are concerned with gravitationally 
stable disks, we need to choose appropriate values for /s 
and /t satisfying the condition given by equation (|13p . 
Assuming the baryonic mass of hyperaccretion disks is 
Mdi s k — lM©, the outermost radius of the disk is about 
70r s for the ADAF-type disk and 30r s for the NDAF-type 
disk. Then, in order that the entire disk is gravitation- 

ally stable, these parameters should be / s f T ' > 1 and 

/e Vt /2 £ 3 for the ADAF-type and NDAF-type disks, 
respectively. We consider the evolution of the disks with 
gravitationally unstable outer regions in § 6. 

3. LINEAR GROWTH OF THE MRI IN HYPERACCRETION 

DISKS 

We evaluate the Reynolds number defined by equa- 
tion (TT]) for hyperaccretion disks using equations (|10[) and 
(fT2")) . The radial profile of the Reynolds number is given 

by 

Re = 2.1 x 10- 7 /^ 3 /- 19/6 /lM 3 5/3 r( 19p - 149 - 12 )/ 6 . 

The Reynolds number is independent of the radius f for 
the ADAF models. For the NDAF models, on the other 
hand, the Reynolds number increases as f increases. 
Then the neutrino viscosity is more effective in the inner 
part of the disks. 

The condition that the growth rate of the MRI takes 
of the order of the angular velocity is Re > 1. For the 
case of /s = fx = 1 (our fiducial model), it requires 
that the magnetic field parameter should be > 10 3 . 
This indicates that the strong magnetic field more than 
10 14 G is necessary for rapid growth of the MRI in hyper- 
accretion disks. When the magnetic field is weaker than 
this critical value, the inner region of the disk could be a 
"dead zone" where the neutrino viscosity suppresses the 
MRI. To show this picture more concretely, in the follow- 
ing sections, we solve the dispersion equation including 
the neutrino transport, and derive the maximum growth 
rate of the MRI in hyperaccretion disks. 

3.1. Dispersion Equation with the Neutrino Viscosity 

MSS07 derive the dispersion equation for the MRI in- 
cluding the effects of the neutrino transport. We apply 
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it to find out the most unstable modes of the MRI in 
hyperaccretion disks. 

The effects of the stratification due to the thermal and 
leptonic gradients can be ignored in the disk. We focus 
on the linear growth rate of the axisymmetric MRI. This 
is because the axisymmetric mode is the fastest grow- 
ing mode of the MRI, and the stability conditions for 
the axisymmetric and nonaxisymmetric MRI are identi- 
cal (MSS07). The toroidal component of the field does 
not affect the linear and nonlinear growth of the MRI if 
the field strength is subthermal (Blaes & Balbus 1994; 
Sano & Stone 2002), so that only the poloidal field is 
considered in our dispersion equation. Ohmic dissipa- 
tion would be negligible in hyperaccretion disks. 

Assuming a uniform vertical field and considering only 
the damping effect by the neutrino viscosity, the disper- 
sion equation in MSS07 can be written as 

7 4 + a 3 7 3 + a 2 7 2 + 017 + a = , (15) 

where 

a 3 = 2vk 2 , a 2 = v 2 k A + 2(k z v A ) 2 + {k z /k) 2 n 2 , 

ai = 2vk 2 {k z v A ) 2 , a = (k z v A ) 4 -4(k z /k) 2 (k z v A ) 2 n 2 . 

Here 7 is the growth rate, k is the epicyclic frequency, 
v is the neutrino viscosity, and k — (k 2 + k 2 ^) 1 ! 2 is 
the wavenumber. Radial and vertical wavenumbers are 
described by k r and k z . For the Keplerian disks, the 
epicyclic frequency is equal to the angular velocity. Note 
that the dispersion equation (fT4|) is characterized by a 
single parameter, which is the Reynolds number Re — 
v A /(vSl). We solve this dispersion equation numerically 
and show the maximum growth rate of the MRI in hy- 
peraccretion disks. 

3.2. The Maximum Growth Rates of MRI 

Figure [1] shows the maximum growth rate of the MRI 
as a function of the disk radius for the cases with different 
magnetic parameters fg = 1, 10, 10 2 and 10 3 . Upper 
and lower panels indicate the results for the ADAF-type 
disk (p = 1.0, q = 0.5) and the NDAF-type disk (p = 
1.5, q — 0.0), respectively. The vertical and horizontal 
axes are normalized by the Keplerian angular velocity 
Q.K and the Schwarzschild radius r s . The parameters /s 
and frp are assumed to be unity in this figure (the fiducial 
model). The critical radius dividing the neutrino-opaque 
and neutrino-transparent regions locates at r cr it ~ 20r s . 
for the ADAF-type disk and r crit « 13r s for the NDAF- 
type disk. 

The maximum growth rate of the MRI in the neutrino- 
opaque region is much smaller than that in the neutrino- 
transparent region if the magnetic parameter fg is 
smaller than 10 3 . That is, the MRI is strongly sup- 
pressed by the neutrino viscosity when the magnetic field 
is weaker than the critical value, B C rit ~ 10 14 G. Turbu- 
lent motions at the nonlinear stage cannot be sustained 
when the growth rate of the MRI is much less than the 
angular velocity fl (Sano & Stone 2002). In contrast, 
MHD turbulence driven by the MRI can grow actively 
in the neutrino-transparent region regardless of the field 
strength, because the viscosity effect can be neglected 
there. 

The characteristics of the MRI in hyperaccretion disks 
depend on the parameters /s and fx- The maximum 



growth rate of the MRI for the cases with the different 
disk temperature and surface density are investigated for 
the two types of disk models in Figures [Hand [31 Normal- 
izations of the vertical and horizontal axes are the same 
as in Figure [TJ The magnetic parameters fg is assumed 
to be unity in these figures. 

For the both disk models, it is found that the maximum 
growth rate of the MRI decreases and the size of the 
dead zone expands with increase of the disk temperature. 
This is because the neutrino depth and viscosity increase 
with the disk temperature [see eqs. ([2]) and ||3J)]. On the 
other hand, the decrease of the surface density yields the 
shrinking of the dead zone and decrease of the maximum 
growth rate. The contraction of the dead zone arises 
from the decrease of the neutrino depth. The reduction 
of the maximum growth rate is caused by the rise of the 
neutrino viscosity with decrease of the surface density. 

These results indicate that the dead zones with the 
size of ~ 10 - 20r s are accompanied with hot and dense 
hyperaccretion disks (/s ~ fx ~ 1). In what follows, we 
develop an evolutionary scenario of hyperaccretion disks 
with the relatively large size of dead zone. 

4. EPISODIC ACCRETION MODEL FOR GRBS 

As is described in the previous section, the MRI is sup- 
pressed by the neutrino viscosity in the neutrino-opaque 
region when the magnetic field is weaker than the crit- 
ical value, £?crit ~ 10 14 G. In contrast, the neutrino- 
transparent region (hereafter active region) is unstable 
for the MRI and its growth rate is of the order of the 
angular velocity. These features of the MRI in hyperac- 
cretion disks have not been considered in previous work 
at all. However, this can be an important physical clue 
to reveal the central engine of GRBs. 

The existence of the dead zone could cause episodic 
hyperaccretion, which can be the source of the short- 
term variability in the prompt emission of GRBs. We 
focus on the hyperaccretion disk formed in the context 
of "collapsar" scenario (Woosley 1993). In this scenario, 
the central ~ 3M of a massive star collapses directly to 
a black hole, while infalling material with large specific 
angular momentum forms an accretion disk of its mass 
~ 1 - 2M Q . 

4.1. Evolutionary Scenario of the Disks with Dead 
Zones 

We consider a relatively large dead zone (r cr it ~ 20r s ) 
formed at the inner part of a collapsar disk. In such 
the case, angular momentum transport in the dead zone 
would be taken by the neutrino viscosity itself, not by 
the turbulent viscosity sustained by the MRI. The size 
of the neutrino viscosity can be described by using the 
a-parameter, 

a„ = ^ = 5.2x 10- 5 f 2 f 2 M- 2 f- 2 ^+ 3 / 2 . (16) 

Thus, the a- value of the dead zone is typically ~ 10~ 4 for 
the fiducial model (/s = fx = 1)- The mass accretion 
rate, M utj onto the central black hole is given by 

M out = 4irrpHv r ~ 7.2 x 10~ 4 (^) hf T M 3 M Q , 

(17) 

where the drift velocity is assumed to be v r = 3v/(2r) 
and M© is the mass accretion rate of IM@ sec -1 . This is 



DEAD ZONE FORMATION AND NONSTEADY HYPERACCRETION IN COLLAPSAR DISKS 



5 



the mass of the baryonic matter extracted from the dead 
zone per unit time. 

On the other hand, in the active region, the turbu- 
lent viscosity driven by the MRI would play an essen- 
tial role in the angular momentum transport. It is ex- 
pected from various nonlinear simulations of the MRI 
that the a-parameter due to the turbulent viscosity could 
be a t ~ or less (Hawley et al. 1995; Stone et al. 
1996; Sano et al. 2004). Therefore, the mass inflow rate 
into the dead zone M- m is about 

M in ~ 7.2 x 10- 2 ( /e/tM 3 M (18) 

This is the mass of the baryonic matter which flows into 
the dead zone per unit time. 

From the equations (|17p and l|18p. it is expected that 
the baryonic matter is accumulated into the dead zone. 
The mass accumulation rate into the dead zone M accu is 
equivalent to the mass inflow rate, because the extracted 
mass from the dead zone is negligible (i.e., M;„ 3> M ut)- 
If the mass accumulation with a rate M accu m M- ln con- 
tinues, the inner dead zone becomes gravitationally un- 
stable at some stage. Then the gravitational torque may 
cause intermittent mass accretion and drive discrete jets, 
which is thought to be the origin of short-term variabil- 
ity in the prompt emission of GRBs (Rees & Meszaros 
1994; Kobayashi et al. 1997). 

The prospective evolutionary scenario of collapsar 
disks with dead zones is depicted schematically in Fig- 
ure |U The dead zone region undergoes the quiescent 
phase and burst phase repeatedly. At the quiescent phase 
(Fig. |4jz), the neutrino viscosity is the only mechanism of 
angular momentum transport. The mass accretion rate 
onto the black hole is relatively small and the mass of 
the dead zone increases during this phase. When the 
gravitational instability set in, the gravitational torque 
would enhance the accretion rate more than two orders 
of magnitude (Fig. \$$>) . The gravitational torque works 
until the dead zone region becomes gravitationally stable 
as a result of the intense mass accretion onto the central 
black hole. 

The evolution of accretion disks with the dead zone 
is also discussed within the paradigm of protoplanetary 
disks (Nakano 1991; Jin 1996; Gammie 1996; Sano & 
Miyama 1999; Sano et al. 2000; Armitage et al. 2001). 
The origin of the dead zone in protoplanetary disks is the 
ohmic dissipation. Because protoplanetary disks are so 
cold that the ionization fraction is very low. The ohmic 
dissipation of the electric current is efficient enough to 
suppress the linear growth of the MRI. In the next sub- 
section, we discuss the typical evolutionary timescales 
and energetics of the episodic accretion in collapsar disks 
referring to the studies for protoplanetary disks. 

4.2. Timescales and Energetics of the Episodic 
Accretion 

We consider an ADAF-type disk (p = 1.0, q = 0.5) 
described by equations (|4|)- piT|) as an example. Then 
the baryonic mass of the inner dead zone is 

Mdcad = J 2^dr ~ 0.15/ s Mf (^) ' M& ■ 

(19) 



Thus, the most of the baryonic matter is stored in the 
outer active region at the early evolutionary stage. Since 
Toomre's Q value of the disk is 

the dead zone region is in the quiescent phase at the be- 
ginning. As time passes, the matter from the outer ac- 
tive region is accumulated near the critical radius. The 
gravitational instability begins to grow when the surface 
density at the critical radius becomes an order of magni- 
tude larger than that of the initial quiescent phase (i.e., 
h >8). 

Once the outer edge of the dead zone becomes gravita- 
tionally unstable, spiral waves would be exited globally 
in the dead zone, and which cause the intense mass ac- 
cretion onto the central black hole. The mass accretion 
rate due to the gravitational torque can be described as 

M g ~ 0.35/ s / T M 3 M , (21) 

where a g = 0.05 is the a-parameter in gravitationally 
unstable disks (e.g., Lodato & Rice 2005). The duration 
of the intense accretion phase should be comparable to 
the viscous timescale in the dead zone (Armitage et al. 
2001); 

^-^°- 93 ^ M3 (c^)(2lfe) 3/2 — ™ 

The gravitational energy released during an intense ac- 
cretion phase is 

E e * ,M gTgC 2 * 2.4xlO- /s Mf (^) 2 (^-f 2 (^) erg , 

(23) 

where rj is the energy conversion efficiency from the rest- 
mass energy to the gravitational one for the case of Kerr 
black holes (Shapiro & Teukolsky 1983). 

After the intense mass accretion phase, the dead zone 
returns to a gravitationally stable state because of the 
decrease of the mass in the dead zone region. However, 
the mass accretion from the outer active region to the 
dead zone continues constantly. The quiescent disk pro- 
gressively evolves to the intense accretion phase again. 
This cycle would be repeated and the explosive energy 
release occurs intermittently until the baryonic matter of 
the outer active region is exhausted. Finally, the dead 
zone disk with a low mass accretion rate would be left 
after the episodic accretion stage (see Fig. Sb) . 

If the episodic accretion is the origin of multiple rel- 
ativistic shells, which are expected to cause the inter- 
nal shocks (Kobayashi et al. 1997), the typical variable 
timescale r var in the prompt emission of GRBs should 
correspond to the duration of an intense accretion phase: 

r var ~ r g ~ 0.93 sec . (24) 

This could be the origin of the observable log-normal fea- 
ture in the short-term variability of the prompt emission 
(Nakar & Piran 2002; Kobayashi et al. 2002; Shen & 
Song 2003). When a few percent of the gravitational 
energy is converted to the radiative one, the typical lu- 
minosity is evaluated as 

L var ~ / A ~ 2.6 x 10 51 (-J-) erg sec' 1 , (25) 
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where / is the conversion factor from the gravitational 
energy to the radiative one. These are almost identical 
to the observed timescale and peak luminosity of variable 
components in the prompt emission (Norris et al. 1996; 
Nakar & Piran 2002). 

The episodic mass accretion is terminated when the 
material in the outer active region is exhausted. Then, 
the total duration of the prompt emission will be deter- 
mined by the mass depletion timescale. Applying the 
mass inflow rate given by equation (jT5J) , we can evaluate 
the total duration and luminosity of the prompt emission 
as 

not ^ Mtot/Mn O(10) sec, (26) 
Lt„t * / * O(10 51 ) erg sec- 1 , (27) 

Ttot 

where M to t is the total mass of the accretion disk and is 
assumed to be 1 - 2 Mq. Thus, our episodic accretion 
model of the collapsar disk can explain many observed 
features of long GRBs quantitatively. 

5. DISCUSSION 

Recent Swift observations have discovered the flaring 
activity in the early X-ray afterglow occurring « 10 - 
10 4 sec after the prompt burst (Burrows et al. 2005; 
Nousek et al. 2006; O'Brien et al. 2006). The flare is 
suggested to be a prolonged activity of the central engine, 
and various models are proposed for explaining its origin 
(Ioka et al. 2005; King et al. 2005; Proga & Zhang 
2006; Perna et al. 2006). How can the flaring activities 
be interpreted in the context of our episodic accretion 
scenario? 

We consider the case that the outermost region of the 
collapsar disk is gravitationally unstable, that is, Q < 1 
around 70r s . When the disk is gravitationally unstable, 
it evolves toward the spirally-undulated structure or frag- 
ments into bound objects. Fragmentation occurs when 
the local cooling time in the disk is shorter than the or- 
bital time ti^ 1 (Gammie 2001; Rice et al. 2003). In this 
case, the local cooling dominates the viscous heating sus- 
tained by the gravitational instability. 

In the outermost region of the collapsar disk, the cool- 
ing of 4 He photodisintegration can be effective (Piro & 
Pfahl 2006; Kawanaka & Mineshige 2007). The cool- 
ing time would be much smaller than the orbital time 
there. Thus, if the outer region of collapsar disk is gravi- 
tationally unstable, it can fragment into multiple bound 
objects. The mass of a bound fragment is estimated as 
rrif ~ 0.2M Q assuming the coalescence of small frag- 
ments into a single body (Piro & Pfahl 2006). 

Considering that the orbital angular momentum of the 
fragment is mainly extracted in the form of gravitational 
radiation, its infalling time is typically 

^-6.5xl03Ml(^) '(^sec , (28) 

(Peters 1964; Piro & Pfahl 2006), which agrees with the 
typical arrival time of the X-ray flares (Burrows et al. 
2005; Nousek et al. 2005; O'Brien et al. 2006). Note 
that this is much longer than the evolutionary timescale 
of hyperaccretion disks discussed in the previous section. 

As the fragment approaches to the central black hole, 
it suffers strong tidal force. Then, the bound fragment 



is disrupted and a rotating disk of its mass ~ 0.2M Q is 
formed around the black hole. The mass of this newly 
formed disk (remnant disk) is comparable to the mass of 
the dead zone region [see eq. (fP9|) ]. Therefore, the entire 
region of the remnant disk is magnetorotationally inac- 
tive if the field strength is weaker than the critical value, 
B crit ks 10 14 G. Unlike the early evolutionary stage of 
the collapsar disk, the episodic accretion can not be ex- 
pected in the remnant disk. This is because the baryonic 
matter does not exist at the outside of the dead zone. 
Therefore, the evolutionary timescale is determined by 
the neutrino viscosity (a v ~ 10 -4 ) and the mass accre- 
tion rate is relatively low, M out ~ 7.2 x lO^Af© [c.f., 
ea(|17p]. The duration of the stationary mass accretion 
would be 

r rem ~m f /M ovlt ~ 0(100) sec, (29) 

which is consistent with the observed duration of the X- 
ray flares (Chincarini et al. 2007). 

6. SUMMARY 

In this paper, we have investigated where the MRI op- 
erates in hyperaccretion disks focusing on the effect of the 
neutrino viscosity in order to reveal the energy releasing 
mechanism of GRBs. Our main findings are summarized 
as follows: 

1. Assuming reasonable disk models and simple neu- 
trino transport, the MRI is suppressed by the neutrino 
viscosity in the inner neutrino-opaque region of hyper- 
accretion disks when the magnetic field is weaker than 
the critical value, -B C rit — 10 14 G. On the other hand, 
in the outer neutrino-transparent region, the interaction 
between the neutrino and baryonic matter can be ne- 
glected, and thus the MRI can drive MHD turbulence 
actively. Most of the disk models show that hyperac- 
cretion disks consist of the inner dead zones and outer 
active regions. 

2. If the MRI is suppressed in the dead zone, the bary- 
onic matter accreted from the outer active region should 
accumulate there. The mass inflow rate from the ac- 
tive region is estimated as M- ln s=s 1O _2 M0. Then, the 
inner part of the disk may become gravitationally un- 
stable at some evolutionary stage. This suggests that 
intense mass accretion onto the central black hole could 
be caused episodically by the gravitational torque. 

3. The typical evolutionary timescale of the intense 
accretion phase corresponds to the viscous timescale of 
the gravitationally unstable dead zone, and is evaluated 
as T g w 1 sec. If the episodic release of the gravita- 
tional energy is the origin of multiple relativistic shells, 
the typical variable timescale and peak luminosity in the 
prompt burst are given by r var w r g ~ 0(1) sec and 
L var ~ 0(1O 51 ) erg sec -1 . These are consistent with 
the observational features of the short-term variability 
in GRBs. 

4. The total duration of the prompt burst should be de- 
termined by the mass depletion timescale of the outer ac- 
tive region, because the episodic mass accretion is termi- 
nated when the material in the active region is exhausted. 
Considering a collapsar disk of mass ~ 1 - 2 Mq, we 
can evaluate the total duration of the prompt burst as 
Ttot — 0(10) sec. In addition, the typical luminosity of 
the prompt burst is obtained as L tot ~ 0(1O 51 ) erg sec -1 
in our episodic accretion model. 
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(a). ADAF-type disk (p=1.0, q=0.5) 
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(b). NDAF-type disk (p=1.5, q=0.) 
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Fig. 1. — The maximum growth rate of the MRI as a function of the disk radius for the cases with different magnetic parameters fs = 1, 
10, 10 2 , and 10 3 . Other arbitrary parameters are fixed to be unity in this figure, that is, /s = /t = 1- The vertical and horizontal axes 
arc normalized by the Keplerian angular velocity flj{ = 2.4 X 10 4 M 3 _1 f -3 / 2 sec -1 and the Schwarzschild radius r a = 8.9 X 10 5 M3 cm, 
respectively. Here the mass of the central black hole is assumed to be Mbh = 3MsMq. Upper panel represents the case of the ADAF-type 
disk (p = 1.0, q = 0.5), and lower panel is the case of the NDAF-typc disk (p = 1.5, q = 0.0). 
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Fig. 2. — The maximum growth rate of the MRI as a function of the disk radius for the ADAF-type disk with different model parameters 
/e and Jt • Normalizations of the vertical and horizontal axes are the same as in Figure [T] The magnetic parameter fg is fixed to be 
unity. Upper panel shows the /^-dependence of the maximum growth rate for the case of = 2. Middle and lower panels are those for 
the cases of /s = 1 and /s = 0.1, respectively. 
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Fig. 3. — The maximum growth rate of the MRI as a function of the disk radius for the NDAF-type disk with different model parameters 
/s and . Normalizations of the vertical and horizontal axes are the same as in Figure [T] The magnetic parameter fs is fixed to be 
unity. Upper panel shows the /^-dependence of the maximum growth rate for the case of fa = 2. Middle and lower panels are those for 
the cases of /e = 1 and /s = 0.1, respectively. 
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Fig. 4. — Schematic pictures of typical three evolutionary stages of a collapsar disk with a dead zone. Panel (a) shows the mass 
accumulation stage. In this stage, the MRI is suppressed by the neutrino viscosity and angular momentum transport is inefficient in the 
inner neutrino-opaque region. On the other hand, in the outer neutrino-transparent region, the MRI can drive MHD turbulence actively 
and the angular momentum is transported efficiently. Then the baryonic matter is accumulated into the dead zone during this stage. 
Panel (b) represents the intense accretion stage. If the baryonic matter is continually accumulated into the dead zone, it should become 
gravitationally unstable at some evolutionary stage and the intense mass accretion is triggered by the gravitational torque. The intense 
accretion would be terminated when the mass of the dead zone becomes small enough to make the region stable for the gravitational 
instability. However, the mass accretion from the outer region continues, and thus the quiescent dead zone approaches guradually to the 
intense accretion stage again. Therefore, the dead zone region alternates between these two stages and the explosive energy release occurs 
intermittently until the baryonic matter of the outer active region is exhausted. Panel (c) depicts the final stage after the episodic mass 
accretion is terminated. The dead zone disk with a low mass accretion rate will be left around the central black hole. Multiple discrete 
jets originated from the episodic mass accretion propagate with a relativistic speed. These multiple relativistic shells might be an origin of 
the short-term variability in the prompt emission of GRBs. 



